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Abstract 

The flavor-singlet component of the g' meson is related to the topological structure of the SU(3) 
gauge field through the chiral anomaly. We perform a 2 -|- 1-flavor lattice QCD calculation and 
demonstrate that the two-point function of a gluonically defined topological charge density after a 
short Yang-Mills gradient flow contains the propagation of the g' meson, by showing that its mass 
in the chiral and continuum limit is consistent with the experimental value. The gluonic correlator 
does not suffer from the contamination of the pion contribution, and the clean signal is obtained 
at significantly lower numerical cost compared to the conventional method with the quark bilinear 
operators. 
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Among other hadrons the t]' meson plays a special role in the study of the topological 
structure of Quantum Chromodynamics (QCD). The rj' meson would be a pseudo Nambu 
Goldstone boson associated with the spontaneous breaking of the axial U{1) symmetry, 
while it acquires a large mass through the chiral anomaly [l|, which relates the divergence 
of the flavor-singlet axial vector current to the topological charge density in QCD. Witten 
and Veneziano ^ estimated the rj' meson mass in the large-W (number of colors) limit 
and showed that its mass squared is proportional to the topological susceptibility of QCD. 

In real QCD with W = 3 and light dynamical quarks, the argument of Witten and 
Veneziano is no longer valid. It is not the r]' meson but the pion that governs the topo¬ 
logical susceptibility. In fact, in our previous lattice QCD simulations where we kept the 
chiral symmetry (nearly) exact j^-^, it was conhrmed that the topological susceptibility 
is proportional to the light sea quark masses as predicted by chiral perturbation theory, 
Xt = ^ denotes the chiral condensate, and rrii denotes the i-th light quark 

mass. In particular, Xt vanishes in the limit of massless up and down quarks, reflecting the 
long-range dynamics of the pion held. 

Then, an interesting question arises: what happens to the t]' meson with W = 3? Since 
the effect of the anomaly is stronger than in the large-A), limit, the 1 ]' meson mass should be 
still generated by the topological huctuation of the gluons. Nevertheless, it must be insen¬ 
sitive to Xt- This implies a non-trivial double-scale structure in the topological excitation 
of gauge held: it creates the p' meson at short distances, while making a tight connection 
to the pion at long distances. The answer to this question could be that the topological 
property of the rj' meson in QCD is hidden inside the pion clouds. 

In this work, by an explicit calculation in 2 -|- 1-havor lattice QCD, we show that the two- 
point function of the topological charge density at short distances gives a mass consistent 


with the experimental value of 
same correlation functions (see 


die 7]' meson mass. Since we have computed Xt using the 
^ for the details), our data clearly show the double-scale 
structure of the topological property of gluons: the long-range physics described by the pion, 
and the short-range (or hrst excited) physics governed by the p' meson. 

Not only is it theoretically interesting, but our work also proposes a practically advan¬ 
tageous method to determine the r]' meson mass. Direct lattice calculation of the r]' meson 
mass has been challenging because one has to include the disconnected quark-line diagram, 
which appears from the Wick contraction of the fermion bilinear operators representing 
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the = 0 flavor-singlet state, as has been done in recent calculations This is 

numerically demanding and statistically very noisy. 

The advantage of the gluonic calculation adopted in this work over the conventional 
fermionic one is two-fold. First, we can avoid the enormous computational cost of stochasti¬ 
cally evaluating the disconnected diagram. The gluonic definition of the topological charge 
density does not require any inversion of the Dirac operator. 

Second, our method avoids the contamination from the pions. In the conventional method 
where one calculates quark connected and disconnected diagram contributions appearing 
from the Wick contraction of quark fields, the pion may propagate in each contribution but 
it cancels between them. A large statistics is required for a sufficient cancellation before 
extracting the 7]' meson propagator. Since the purely gluonic definition of the topological 
charge density operator 


q{x) = 


1 


.TrF,rF,r(x), 


( 1 ) 


32^2 "/..pc. 

where denotes the held strength tensor of the gluon held defined through the so-called 
clover term consisting of four plaquettes, does not directly couple to the pions, its correlator 
is free from the pion background. 

Note here that the sum of ([T]) over the lattice volume gives the global topological charge 
up to discretization effects, Q = + 0{a^). However, it is well-known that the 

O(a^) contribution is large with currently available lattice spacings. In order to reduce this 
statistical noise, we modify the link variables using the Yang-Mills (YM) gradient flow 12|. 
At a flow time t, it amounts to smearing the gauge fields in a range of the length \/M. It 
is shown that the topological char ge Q defined through ([T]) converges to an integer value 


13 


14| . This smearing procedure eliminates short-distance 


at a sufficiently large flow time 
fluctuations and also suppresses the noise at longer distances.^ 

For our purpose of extracting the r( meson mass, the YM gradient flow time has to 
be short in order not to distort the correlation of the r/' propagation. Assuming a simple 
Gaussian form of the smearing effect, Bruno et al. estimated the size of distortion of 


^ A similar method was applied in a quenched study to extract the “pseudoscalar glueball mass” [l^. There 
are other viable definitions of the smearing as used in previous works to probe the topological structure 
of the QCD vacuum 17-2^. 
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the correlator as 


A{q{x)q{y)) y^)2j7v(8t)^ 


( 2 ) 


{q{x)q{y)) " 2y/n\x - y]"^' 

In our analysis below, we use the reference flow time around \/^ = 0.2 fm for the £t range 
\x — y\ > 0.6 fm, which makes the above correction less than 1% for rurf' — 1 GeV. 


We employ the Symanzik 
gauge ensemble generations 


gauge action and the Mobius domain-wall fermion action for 


21 


23l |. We apply three steps of stout smearing of the gauge 


links before inserting it in the Dirac operator. Our main runs of 2 -|- 1-flavor lattice QCD 
simulations are carried out on two different lattice volumes xT = 32^ x 64 and 48^ x 96, 
for which we set /3 = A.17 and 4.35, respectively. The inverse lattice spacing 1/a is estimated 
to be 2.4 GeV (for (3 = 4.17) and 3.6 GeV (for {3 = 4.35), using the input = 0.1465 fm 
where the reference YM gradient flow time to is dehned by t^{E)\t=tQ = 0.3 12|, with 
the energy density E of the gluon held. Our two lattices share a similar physical size L ~ 2.6 
fm. For the quark mass, we use two values of the strange quark mass rrts around its physical 
point, and 3-4 values of the up and down quark mass rriud for each m^. The lightest pion 
mass is around 230 MeV with our smallest value of arriud = 0.0035 aX (3 = 4.17. In order 
to check the systematics due to hnite volume sizes and lattice spacings, we also perform 
simulations on a larger lattice 48^ x 96 (at (3 = A.17 and ~ 230 MeV), and a finer lattice 
64^ X 128 (at (3 = 4.47 [1/a ~ 4.5 GeV] and tUt, ~ 285 MeV). For each ensemble, 500-1,000 
gauge configurations are sampled from 10,000 molecular dynamics (MD) time. The residual 


2d by 


mass in the Mobius domain-wall fermion formalism is kept smaller than ~ 0.5 MeV 
choosing = 12 at /? = 4.17 and = 8 at /3 = 4.35 (and 4.47). 

On each generated configuration, we perform 500-1,000 steps of the YM gradient how 
(using the conventional Wilson gauge action) with a step-size a^Af =0.01. At every 20- 
30 steps, we store q{x) and measure its correlator using the Fast Fourier Transform (FFT) 
technique. As reported in 6j, the how time history of the gluonic dehnition of the topological 
charge Q shows a good convergence to discrete values near integers. 

We hnd that the two-point function {q{x)q{y)) at our target distance \x — y\ ~ 0.7 fm 
always shows a shorter autocorrelation time than 10 MD time, while its global average, 
Q = J2x has 0(100) or hig her MD time at /? = 4.35. This is a good evidence that the rj' 


meson physics is decoupled 


26(1 from the global topological charge. In the following analysis. 
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FIG. 1: The correlator —{q{x)q{y)) at the flow times a^t = 0.2 (circles), 0.4 (triangles) and 1.0 
(crosses). Data at /3 = 4.17, arUud = 0.007 and anig = 0.040 are presented. 


we estimate the statistical error by the jackknife method after binning the data in 140-200 
MD time. 

Figure [H shows the topological charge density correlator C{\x — y\) = —{q{x)q{y)) at the 
flow times a^t = 0.2, 0.4 and 1.0. Using the FFT, data points are averaged over all possible 
combinations of x and y giving the same r = \x — y\io improve the signal. As the flow-time 
increases, the statistical fluctuation of the correlator becomes milder, while the region at 
small \x — y\ is distorted. We therefore need to hnd a region of t where the correlator has 
sufficiently small noises to hnd the signal while it is not distorted by the smearing of the 
YM gradient how. 

The data for C{\x — y\) show a cleaner signal than the conventional zero-momentum 
projection, '^{q{x)q{y)), because of the average in all four-dimensional directions. The rj' 

X 

meson mass is extracted by htting the data to the function of a single boson propagation: 


/(r,m^0 = A 


Ki(mn'r) 


( 3 ) 


where r = \x — y\, Ki is the modihed Bessel function and A is an unknown constant, which 
depends on the how time t. 

In order to determine the htting range, we dehne a local “ehective mass” mefr(r) by 
numerically solving /(r -|- Ar, mefr(7’))//(r, mefr(r)) = C{r + Ar)/C{r). We set the interval 
to be Ar = 1/2, and the data of C{r) in the range [r, r -|- Ar] are averaged. As shown 
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FIG. 2: The effective mass meff(r) for the data with several flow times. Results for the /3 = 4.17 
= 0.007 and arus = 0.040 ensemble are shown. 
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FIG. 3: The flow-time dependence of the t]' meson mass. The data at various sea quark masses 
and /3 values are shown, as specified in the legend. The filled symbols represent our data taken for 
the central values. 

in Fig. [21 a reasonable plateau is found for meff(r) around r ~ 8-12 (> 0.6 fm) at f = 1 
~ 0.2 fm). 

Figure [3] shows the t]' meson mass obtained by htting in the range [8,12] as a function of 
\/8f- The data around \/8f ~ 0.2 fm are stable. At larger smearing lengths y/St > 0.3 fm, 
we observe a large distortion of the data. We take the data at y/^ = 0.2-0.25 fm (hlled 


6 

















































FIG. 4: The extracted rj' meson mass from each gauge ensemble. The three fit lines (representing 
the same linear fit function at three values of a) are shown for higher rus’s at each /3. 


symbols in Fig. [3]) for our results. 

We plot all the results in Fig. 0] as a function of the square of the pion mass m^. We 
hud that the dependence on the quark masses, as well as on V and a is mild. We therefore 
perform a global £t of our data to a linear function + CaO^ + Cudimi — + 

Cs[{2mj^ — ml.) — where Ca, Cud, and Cg are taken as free 

parameters, and mP^^ denotes the experimental value of the pion/kaon mass. As shown by 
the lines (which are shown for higher mg only) in Fig. 01 we hnd that our linear function hts 
the lattice data reasonably with x^/(degrees of freedom) ~ 1.6. 

Because of possible bias in the topological charge sampling, the rj' correlator may not 
decay exponentially but become a constant at long distances |27| . In each topological sector 


Q, it is predicted as 


28| 


'Q 


{(l{x)q{y)r ~ ^ 1 17 “ I |x - y\. 


(4) 


The typical size of the constant term Xt/V is 100-1000 times smaller than \{q{x)q{y))\ in 
our htting range, which is consistent with the fact that this observable shows no strong 
correlation to the global topological charge Q. 

We also estimate that the hnite volume effect on our observable is negligible, since the 
7]' meson propagator implies exp(—mr,'T) ~ 3 x 10“® even at our lightest mud- Our data on 
a larger lattice 48^ x 96, which are statistically consistent with those on the smaller lattice. 
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FIG. 5: Comparison of the fit curves at different input 9. 


support this assumption. 

Since rriud and rUg are non-degenerate, the flavor-singlet channel does not correspond to 
the mass eigenstate and a mixing with the flavor non-singlet channel is expected, which 
is the rj-T]' mixing. We estimate the size of the p-meson contribution to the flavor-singlet 
channel by modifying the fit function as 


A 


A r 


Kiiniri'r) rurjiKi{mri'r) + tan^ 6 




J- L 


(5) 


where 9 is the mixing angle, and denotes the r] meson mass. Note that the factors 
and come from the effect of the YM gradient flow smearing. As a phenomenological 
estimate for m,,, we use ^{mud + “^ms) /(Sm^) where is the funphysical) mass of the 
connected pseudoscalar correlator of strange and anti-strange quarks 29|. We compare the 
fit curves with different fixed values of 9 as shown in Fig. [5l As 0 increases the quality of the 
fit becomes worse, especially at long distances. Thus, our data suggest a small mixing angle 
\9\ < 10°, and this tendency is seen at all our simulated parameters, (the sign of 9 is not 


relevant for this conclusion.) The sma^ 


1 mixing angle is consistent with the estimates from 


the quark model, 9 = —25°-10° 30|. In our simulations, 9 is likely to be smaller than 


these phenomenological values since our simulated ruud is closer to m^. We confirm that the 
fit results with the fixed value of 9 = 10° are consistent with those with 9 = 0 within the 
statistical errors, although there exists a tendency that the m,,/ meson mass becomes higher 
for 9 0. We take this ~ +5% deviation as the systematic error from the r] meson mixing. 


8 



















Finally, we examine the systematics in the chiral and continuum extrapolation. Since 
both of the and dependences are mild, even if we totally ignore these dependences, 
the (constant) £t works well, giving a 8% different value of from the original linear 

£t, which is within the statistical error. We take this ±8% as the possible systematic error 
in the extrapolations. 

Our final result at the physical point is 


= 1019(119) MeV, 


( 6 ) 


which is consistent with the experimental value mrj> = 957.78(6) MeV 30|. Here the first 
error is statistical and the second is the systematic error from the mixing with the rj meson 
and the chiral and continuum extrapolations (added in quadrature). From the same set of 
correlators we obtain the topological susceptibility Xt, which will be presented elsewhere. 
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